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Abstract 

The values of quark condensate and tensor susceptibility of QCD vac- 
uum are calculated with different number of flavors, ft is shown that, 
the values of tensor susceptibility are obviously dependent on the flavor 
numbers while the flavor-number dependence of quark condensate is un- 
determined. 
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The quark condensate and tensor susceptibility are two quantities of reflect- 
ing the nonperturbative aspects of QCD vacuum. The recent calculations show 
that, there exists a great discrepancy between our values of the QCD vacuum 
tensor susceptibilities and those obtained from QCD sum rules and from chiral 
constituent quark model J^. Within the framework of chiral perturbation the- 
ory [2], Moussallam argued that the quark condensate is strongly dependent on 
the flavor- number Nj via the calculation of the chiral coupling-constant Lq 
This strong Nf dependence cannot be derived in our calculations below, while 
on the other hand, we find the strong Nf dependence of the values of tensor 
susceptibilities. Our calculations are performed within the framework of an ef- 
fective QCD model, the global color symmetry model (GCM) 0], which has 
been well developed in the last twenty years 0. 

The GCM partition function for massless quarks in Euclidean space is writ- 
ten as 

2gcm = / r'qr»gr'Ae-^'=^<="[«'«"^i (i) 



with the GCM action 

1 



SGCM[q,q,A]^ J dx[q{x){^^ig4)q{x) + J dxdy^A';,{x)[D';:l{x ~ y)]-'Al{y). 

(2) 

Here D'^'l{x—y) is the gluon 2-point Green's function, and a Feynman-like gauge 
^"fiti^ ^ y) ~ ^tiu^°'''D{x — y) is used. The important dynamical characteristics 
of local color symmetry are included in D{x — y). The phenomenological gluon 
propagator is treated as the model input parameter and is adjusted to reproduce 
the pion decay constant. 

Applying the standard bosonization procedure, the resulting expression for 
the GCM partition function in terms of the bilocal field is given by 

^GCM = j VB' exp {-S[B']) , (3) 

where the action is 

S[B]^-T.Ln[G ] + J dxdy ^^.^(^ _ ■ (4) 
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The quark inverse Green's function G ^ is defined as 

G'\x,y) = ^5{x-y) + K'B\x,y), (5) 

where the quantity = \K''C'>F'' arises from Fierz reordering of the current- 
current interaction 



^%Kk = {l,^)jk{l,^)H. (6) 



Here the Dirac matrix is K"- = {lojiln, r^ln, r^luls), the color matrix is 
= (|/c, ^A"), and the flavor matrices for Np = 2 flavors and Np = 3 
flavors are respectively F^^^^ = {:^If, 7^ 'Tf) and F^^^^ = {-^If, -^>'f)- 
The saddle point of the action is obtained by minimizing the bilocal action: 
= 0, which gives 



SS[B] 
SB 



B'oix) = g'D{xMGo{x)A'] . (7) 

We define S(p) = A^B^{p) = ifi[A{p'^) - 1] + B{p'^), where the self-energy 
functions A and B satisfy the Dyson-Schwinger equations, 



[A{p )-l]p = 3 y 9 D{{p-q) 



(27r)4 ^ VV/- ^ 52^2(^2)^52(^2) ' 

4„ Uf„2\ 



(8) 



For the convenience of computation, the Dyson-Schwinger equations are often 
rewritten as 

[A{s)-l]s = -^j^ Als'l ^^n^X^^^(^.^0 ^^,^.(,^,)V^^.(;) dx(10) 



where s = p"^, and g^D{s, s') = g^D{s + s' — 2\/ss' cosx)- 
The quark Green's function at Bq is given by 

-i^Aip^) + B(j)'^) 



Go{x,y) = Go{x-y) = J 



02^2(p2)+52(p2) 

and the equation used to calculate f.^ is 



,ip-{x-y) 



167r2 Jo [sA^{s) + B2(s)]2 [ 5^2(5) + 52(g) 



(12) 



(13) 
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In Eq. all those terms involving the derivatives of A{s) and B{s) with 

respect to s are neglected. Let A{s) = 1, B{s) = M, with M the mass of 
the constituent quark, the integrand of Eq. behaves like 2M^/s when s 
approaches infinity, which reproduce the result of Ref. 6 strictly. 
The two quantities to be evaluated are quark condensate 

/- y ^ r ^ Bis) 

^''^^-I^Jo '''' sAHs) + BHs) ^ 

and the quantity 11^^ (0) defined as 



I^n,(0).-^/ .d. 



(15) 



3A'^{s) + B^{s) 

which relates the quark condensate to the tensor susceptibility x via the defini- 
tion 

We choose three different gluon propagators as follows, the ultraviolet behav- 
ior of them are different from that in QCD |7!IS^. These models are numbered 
as model 1: 

g^D{s)^37T'^e-^ , (17) 



model 2: 



and model 3: 



,^D(.) = 3.^-e- + -^-^-^^, (18) 

g'Dis) = ^^^d^^ . (19) 

Here d = 12/(33 - 2iV/) and A = 200 MeV. 

The calculation results for model 1 are displayed in Tabled and those for 
model 2 and model 3 are respectively given in Tables |5] and |3| It is indicated 
that, the value of in 2-flavor QCD is different from the one in 3-flavor QCD [^. 
However, this difference is so small that almost we need not make any modifi- 
cations to the calculation results. In the model 3 for example, when N f — 2, 
if we take A = lO^^ QcV*, A ^ 1.355 GeV, is fitted to 88 MeV, the values 
of ~{qq)'^l'^ and -nx(0)/12 are respectively 245 MeV and 1.51 x 10^ MeV^ 



4 



which are very close to the corresponding data 243 MeV and 1.47 x 10^ MeV^ 
for ^ 87 MeV. So we take U = 87 MeV in all cases. 

Let's analyse the results listed in Table^for model 1. If A is independent on 
A^/, say, A = 0.200 GeV^ then the parameter A = 1.509 GeV for A^/ = 2 should 
be replaced by A = 1.309 GeV for N f = 3 to fit the pion decay constant, which 
means the parameter A is dependent on 7V^. Similarly, if A is independent 
on A^/, then A must be dependent on (the calculations for this case are 
not given here for clarity). So, in any case, there is at least one parameter 
which is dependent on N f. Obviously, if A is independent on A^/, the values of 
quark condensate —{qq)^^'^ are universally decreased by more than 10% when 
Nf changes from 2 to 3. However, the Nj dependence of A in GCM is actually 
unknown. On the other hand, when A moves from 0.200 GeV^ to 0.500 GeV^ 
the quark condensate -{qqY^^ for Nf = 2 ranges from 207 MeV to 243 MeV. 
If A < 0.200 GeV^ or A > 0.500 GeV^, the corresponding values of quark 
condensate will leave this range. The values of quark condensate for Nf — 3 
ranges from 180 MeV to 216 MeV, which partly coincides with the one for 
Nf — 2. Consequently the Nf dependence of quark condensate is unsure. In 
contrast to this, the Nf dependence of tensor susceptibilities is evident. The 
values of-nx(0)/12 for Nf ^2 ranges from 1.71 x 10^ MeV^ to 2.07x10^ MeV^ 
while that for Nf decrease to (1.25 1.46) x 10^ MeV^ 

It is easy to see that, the similar peculiarities for the values of ~{qq)^^^ and 
— n^(0)/12 can also be found both in Tables |21 and |3| 

In conclusion, for three models of gluon propagators, the values of —{qq)^^^ 
and — n^(0)/12 are obtained by adjusting the parameters A and A to fit the 
pion decay constant. Because the accurate flavor-number dependence of the 
gluon propagators (i.e. the Nf dependence of A and A) is unclear, the flavor- 
number dependence of two quark condensate is undetermined. Compared with 
the results for A'^ — 2 flavors, the values of — n^(0)/12 for A'^ = 3 flavors 
are universally lowered. From Eq. 1161 the definition of tensor susceptibility, for 
an average value of quark condensate, a smaller value of — 11^(0) /12 implies a 
smaller value of tensor susceptibility. The tensor susceptibilities of QCD vacuum 
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are relevant with the flavor numbers. 
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Table 1: {qq) and U^iO) are calculated for model 1: g'^D{s) = Stt^^c-^. The 
parameters A and A are adjusted to fit = 87 MeV. 







Nf = 2 








iVF = 3 




A 

[GeV^] 


A 

[GeV] 


[MeV] 


-n^(0)/i2 

[xlO^ MeV^] 


A 

[GeV^] 


A 

[GeV] 


[MeV] 


-nx(o)/i2 

[xl03 MeV^] 


0.200 


1.509 


207 


1.71 


0.200 


1.309 


180 


1.25 


0.300 


1.604 


221 


1.87 


0.300 


1.439 


194 


1.35 


0.400 


1.709 


233 


1.98 


0.400 


1.568 


206 


1.42 


0.500 


1.816 


243 


2.07 


0.500 


1.690 


216 


1.46 


Table 2: Same as Table [T] for model 2: g'^D{s) = ^■n'^^e^^ ^ 
d ^ 12/(33 - 2Nf), A ^ 200 MeV. 


sln[s/A^+e] ' 








Nf = 2 








NF = i 




A 

[GeV'] 


A 

[GeV] 


[MeV] 


-nx(o)/i2 

[xlO^ MeV^] 


A 

[GeV^] 


A 

[GeV] 


[MeV] 


-nx(o)/i2 

[xl03 MeV^] 


0.100 


1.453 


216 


1.33 


0.100 


1.213 


188 


0.925 


0.200 


1.509 


232 


1.52 


0.200 


1.306 


205 


1.07 


0.300 


1.593 


245 


1.66 


0.300 


1.420 


218 


1.16 


0.400 


1.687 


256 


1.76 


0.400 


1.536 


228 


1.22 


Table 3: 


Same as Tablen]for 


model 3: g'^D{s) 


= 47r2d^,d= 12/(33- 2iV;). 








Nf = 2 








NF = i 




A 

[GeV^] 


A 

[GeV] 


[MeV] 


-n^(0)/i2 

[xlO^ MeV^] 


A 

[GeV^] 


A 

[GeV] 


[MeV] 


-n^(0)/i2 

[xl03 MeV^] 


10-2 


1.345 


243 


1.47 


10-2 


1.157 


207 


1.03 


10-3 


1.109 


223 


1.34 


10-3 


0.930 


187 


0.914 


lo-'^ 


0.955 


210 


1.29 


lo-'^ 


0.791 


175 


0.864 


10-5 


0.847 


201 


1.28 


10-5 


0.695 


167 


0.850 
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